arXiv: 1501.04619vl [astro-ph.HE] 19 Jan 2015 


Maximum mass of magnetic white dwarfs 


D. Manreza Paref] 

Facultad de Fisica, Universidad de la Habana, San Ldzaro y L, Vedado La Habana, 10400, Cuba 

A. Perez Martinet!] 

Instituto de Ciberndtica Matemdtica y Fisica (ICIMAF) Calle E esq 15 No. 309 Vedado 

La Habana, 10400, Cuba 

J. E. HorvatlfO 

Instituto de Astronomia, Geofisica e Ciencias Atmosfericas USP, Rua do Matao 1226, 05508-900 Sao Paulo SP, Brazil 

(Dated: January 20, 2015) 

We revisit in this work the problem of the maximum masses of magnetized White Dwarfs (WD). The impact 
of a strong magnetic field onto the structure equations is addressed. The pressures become anisotropic due to 
the presence of the magnetic field and split into a parallel and perpendicular components. We first construct 
stable solutions of TOV equations for the parallel pressures, and found that physical solutions vanish for the 
perpendicular pressure when B > 10 13 G. This fact establishes an upper bound for a magnetic field and the 
stability of the configurations in the (quasi) spherical approximation. Our findings also indicate that it is not 
possible to obtain stable magnetized WD with super Chandrasekhar masses because the values of the magnetic 
field needed for them are higher than this bound. To proceed into the anisotropic regime, we derived structure 
equations appropriated for a cylindrical metric with anisotropic pressures £T). From the solutions of the structure 
equations in cylindrical symmetry we have confirmed the same bound for B ~ 10 13 G, since beyond this value 
no physical solutions are possible. Our tentative conclusion is that massive WD, with masses well beyond the 
Chandrasekhar limit do not constitute stable solutions and should not exist. 
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I. INTRODUCTION 


Motivated by the observation of thermonuclear supernovae that seem to require exploding white dwarf masses above the 
celebrated Chandrasekhar limit 0, a series of papers by Mukhopadhyay and collaborators J3] (4J explored the magnetized 
version of the stellar structure and argued for a substantial increase of the maximum possible mass for large values of the 
magnetic field B, which quantizes the electronic energy levels. A great deal of interest followed this suggestion, which has been 
addressed in a number of works. The main criticisms include an inconsistency with the virial theorem 0 for the achievement 
of large maximum mass values and similar basic facts. However, the authors have responded to these criticisms, and the issue of 
the existence of these compact stars is still an open and important question. 

On the theoretical point of view, the construction of fully consistent equilibrium solutions in the magnetized regime is still 
lacking, although hints for their existence and stability have been pointed out 0. The spatial distribution of the magnetic field 
seems to be an important ingredient for this issue, while the behavior of matter under extreme conditions leads to a consideration 
of the equation of state in the Landau regime for the electron energy levels, which may change the effective description in 
terms of polytropic indexes and related quantities. Therefore, a step forward towards the solution of this problem would be 
the investigation of the stability of stellar models for highly magnetized matter and identify the threshold values of B for the 
disappearance of stable solutions. We perform in this Letter such an analysis within a definite relativistic framework and show 
that, at least within these simplified models, the magnetic field admitted on theoretical grounds can not exceed 1.5 x 10 13 G. 
Moreover, we confirm that the maximum masses do not grow beyond the Chandrasekhar value when the magnetic pressure 
is properly introduced via the stress-energy tensor. A brief discussion of the theoretical situation of putative high-mass white 
dwarfs closes this work. 
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II. MODELS OF MAGNETIZED WHITE DWARFS 


The scalar virial theorem J6| has been generally employed to estimate that the maximum magnetic field that a WD can sustain 
with a mass M = IAMq and a radius R = O.OO5f?0 is around B mm ~ 10 13 G. This values strongly suggest that a realistic 
model of a magnetized WD should feature quantized energy levels for the electrons, as in many attempts to construct models 
that describe the microphysics of WDs as a magnetized fermion systems |7). 

In the approximation in which the magnetic field is constant and matter is allowed to settle in it, a breaking of the spherical 
symmetry of the star is apparent. This is not very relevant for low magnetic fields, but since we want to reach the extreme 
anisotropic regime, we have chosen to work in cylindrical coordinates in which the polar and equatorial radii differ and the 
deviation of spherical symmetry is naturally accounted for. An additional advantage of this procedure is that the construction of 
an anisotropic energy momentum tensor for the magnetized matter is very well-defined and straightforward. 

In 01 we have developed a first attempt to investigate this problem using a general cylindrical symmetric metric, whit co¬ 
ordinates ( t , r, <f>, z). We have followed there the procedures of ID to solve Einstein equations for an axisymmetric model of a 
WD to take into account the anisotropy induced by the magnetic field. A constant magnetic field in (say) z-direction defines two 
main directions in space, parallel and perpendicular to the magnetic field. The main approximation performed in that paper is to 
assume that all the functions and variables depend only on the radial coordinates (r) and not on z and </>, so that we can solve 
the dependence on the equatorial direction of the WD. However, this simple model could be useful to obtain information of the 
effects of the magnetic held in terms of the shape (oblateness) of the WD and yield upper limits for the values of the magnetic 
held that this objects can sustain. 

The present paper builds on the ID results, and applies the same procedure to study the structure equation of a magnetized 
WD with the aim to conhrm or dismiss the recently claims of super-Chandrasekhar masses for magnetized WD 0. 


III. MAGNETIZED WHITE DWARFS 


The thermodynamical properties of matter in a magnetic held are obtained starting from the thermodynamical potential at 
zero temperature 
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where fi e is the electron chemical potential, l ma x = I[ z \ denotes the integer part of z, ai = 2 — Sio is the spin 

degeneracy of the (-Landau level, the Fermi momentum is pp = \JHe' 2 - £ 2 and the rest energy is given by 
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The particle number density and the magnetization are 
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where B° = m 2 /\e\ = 4 x 10 13 G is the critical magnetic held and C e = -^1/^21-^ + m 2 . 

The energy density and the pressures parallel and perpendicular to the magnetic held can be written as 


e = n e + n e N e + Nm N — , 
'Pw = 

V± = — fl e — BMe , 
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where Nuin ^ is the mass density term, N is the number of nucleons, m.v is the mass of nucleons, Z is the atomic number and 
A baryon number. We assume that the white dwarfs are predominantly composed by 12 C and le O with A/Z = 2. 
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The EoS including matter and field (Pf = E B = — P B = 

E = 

P\\ = 

P± = 


tt-) contributions have the following form 
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FIG. 1: Equations of state for magnetized electron gas. Notice the differences of the pressures when the magnetic field increase, the cases of 
B = 0 and B = 10 12 almost can not be distinguished one of the other. 

In Figure [T] we show the EOS of the magnetized gas as derived from the above expressions. It can be noticed that for low 
values of the magnetic field (B -C /j';), the pressures do not differ much among them and from the non-magnetic case B = 0. 
However, when B ~ B)';, the difference of the pressures becomes quite large. 

To quantify the anisotropy we have defined the splitting coefficient as 


A = 


P(B -> 0)' 


(9) 


We will use as a criterion that the border separating the isotropic and the anisotropic regions is just A ~ 0(1). We have 
solved the equation A(/z e , B) = 1 in order to distinguish the anisotropic region from the isotropic one. In Figure [5] we show 
the plane of matter densities as a function of the magnetic field, the region above the dashed curve falls in the isotropic regime 
(A < 1) and the grey region below the curve is in the anisotropic regime (A > 1). The vertical lines are the solutions of the 
densities for a constant magnetic field value, for B < 10 12 G, all the points lay in the isotropic region while for B > 10 13 G 
there is a considerable number of points in the anisotropic region. 

In our numerical computations we will first use magnetic field values well within the isotropic region B = 10 12 G and also 
in the anisotropic region B = 10 13 G to compare the effects on the star structure. In next section it will be shown the impact on 
the stellar structure of this density-dependent field anisotropy. 
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FIG. 2: Splitting of the pressures whit respect to the magnetic field. The dashed line represent the solution of the equation A = 1 and 
the vertical lines are solutions of the densities for a constant magnetic field value. 


IV. TOV EQUATIONS FOR MAGNETIZED WHITE DWARFS 


In order to set up the problem introduced by the magnetized matter EoS in the study of the structure of WDs, we will analyze 
the usual case first, assuming spherical symmetry and solving the resulting TOV equations (9)- To find the static structure of a 
relativistic spherical star we have to solve the well-known TOV equations. 
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with the boundary conditions P(R) = 0, M(0) = 0 and the EoS E —> f(P). 

The Mass-Radius curves obtained for magnetic fields B = 10 11 ,10 12 G (that is, within the regime in which A < 1) are in 
agreement with the ones obtained, for example, in ifTOl . Therefore, we confirm that quantization of the electronic levels for fields 
in this regime can not increase the maximum mass of a WD sequence. However, if the magnetic field pressure terms oc B 2 /8tt 
are omitted, higher values can be achieved. We believe that this unjustified omission is a significative part of the discussion on 
super- Chandrasekhar masses. 

When we want to explore the anisotropic region, however, we find that the anisotropy sets in for progressively large regions 
of the WD when the value of the magnetic field is increased. Around B = 10 13 G most of the star feels the anisotropy of the 
pressures and just the inner regions remain practically isotropic. This justifies the use of anisotropic solutions for the highest 
fields. 


V. ANISOTROPIC STRUCTURE EQUATIONS 


In order to improve the structure equations in presence of anisotropic pressures we consider in this section an axisymmetric 
geometry which is more adequate to treat magnetized fermion system. We follow the same procedure as in m. The cylindrically 
symmetric metric reads 


ds 2 = — 


e 2 *dt 2 + e 2A dr 2 + r 2 # 2 + e 2 *dz 2 


( 12 ) 







where <t\ A, Cl, and 'I' are functions of r only which as mention before is a main approximation [Q]. 
With the energy momentum tensor for magnetized mater given by 0 
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where E, P\\ and P± are given by the EoS ([6]), dTJi and (|8j) respectively. 

From the Einstein field equations in naturalunits and using the energy momentum conservation we obtain the 

following four differential equations: 


P’ ± = -<i> , (E + P ± )-^(P ± -P ll ), (14a) 

<T>' 

47re 2A (P + P||+2Pj_) = <f>" + <f>'(4'' + $' - A') + —, (14b) 

47re 2A (P + P||-2Pj_) = -tf" _ _ A')-, (14c) 

47re 2A (P|| — E) = -(4'' + $' — A'). (14d) 

This, together with the EoS E —> f(P±), P\\ —> f{E) is a system of differential equations in the variables 

Pl, P||, E, <b, A, vp, (15) 


We consider P±(R±) = 0 which determines the radii of the star, in the equatorial (perpendicular) direction. The solution of 
the system of equations ( fl4| are shown in Figure [3] 
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FIG. 3: Mass per unit parallel length (M/R \\) in solar mases, versus perpendicular radius. We have plotted curves for the following magnetic 
fields values: B = 0 G (isotropic case), B = 10 11 G, 5 x 10 11 G, 10 12 G 5 x 10 12 G and B = 10 13 G. This last value of the magnetic field 
is the maximum value that we found stable configurations. 

We have plotted in Figure [3] the mass per unit length versus perpendicular radius, (we can not calculate a standard Mass- 
Radius relation because we did not retain a dependence of our magnitudes on the ^-direction for this choice of the metric). 
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When the magnetic field increases, the perpendicular radius and the mass per unit length of the star also increase. We have also 
found that there is a maximum field (B ~ 10 13 G) beyond which the metric coefficients exhibit a divergent behavior, this value 
of the magnetic field coincides with the value when the splitting of the pressures given by the parameter A is greater than 1, for 
B = 10 13 G. This result supports our interpretation that beyond this value of the magnetic field there are no stable solutions of 
the system, and point towards the end of the theoretical stellar sequences constructed from our assumptions. 


VI. CONCLUSIONS 

We have revisited the role of anisotropic pressures in the description of the structure of a WD. Our findings show that when 
the splitting coefficient A > 1, the differences in the pressures can not be neglected and a different approach must be used 
to study the structure of the star. An axisymmetric geometry is more suitable than a spherical one for the solution of Einstein 
equation using a cylindrical symmetric metric. Our choice of the metric in the latter conditions is probably the simplest among 
all possible cylindrical metrics. 

Taking into account the pressure anisotropy due to a magnetic field yields a critical field B c ~ 10 13 G for magnetized WD, 
beyond which there are no stable equilibrium configurations. This bound for the value of the magnetic field is close (but slightly 
lower) to the one obtained based on the scalar virial theorem. 

Although in our model we can not compute the total mass due to the assumption that all the variables depend only on the 
perpendicular (equatorial) radius and not of the 2 -direction fT[, the study is useful to confirm the existence of a maximum 
magnetic field for which the star may undergo an anisotropic collapse due to a magnetic instability. This point helps to clarify 
the claim of super- Chandrasekhar masses for magnetized WD 00 and rules out that the magnetic field could be the reason of 
the existence of this kind of objects. By the way, the recent paper ED makes use of extremely high magnetic fields, well above 
the Schwinger value, and clearly beyond the virial estimate; and also imposes a T = 4/3 polytrope as a model for the matter. 
This is at odds with previous claims 0 by the same authors that a T = 2 results from Landau quantization and hence it is clear 
that the latter does not stiffen the equation of state needed to achieve the super-Chandrasekhar mass values. 
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